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We construct a novel scattering theory to investigate magnetoelectricaUy induced spin polariza- 
tions. Local spin polarizations generated by electric currents passing through a spin-orbit coupled 
mesoscopic system are measured by an external probe. The electrochemical and spin-dependent 
chemical potentials on the probe are controllable and tuned to values ensuring that neither charge 
nor spin current flow between the system and the probe, on time-average. For the relevant case of 
a single-channel probe, we find that the resulting potentials are exactly independent of the trans- 
parency of the contact between the probe and the system. Assuming that spin relaxation processes 
are absent in the probe, we therefore identify the local spin-dependent potentials in the sample at 
the probe position, and hence the local current-induced spin polarization, with the spin-dependent 
potentials in the probe itself. The statistics of these local chemical potentials is calculated within 
random matrix theory. While they vanish on spatial and mesoscopic average, they exhibit large 
fluctuations, and we show that single systems typically have spin polarizations exceeding all known 
current-induced spin polarizations by a parametrically large factor. Our theory allows to calculate 
quantum correlations between spin polarizations inside the sample and spin currents flowing out 
of it. We show that these large polarizations correlate only weakly with spin currents in external 
leads, and that only a fraction of them can be converted into a spin current in the linear regime 
of transport, which is consistent with the mesoscopic universality of spin conductance fluctuations. 
We numerically confirm the theory. 

PACS numbers: 72.25.Dc, 73.23.-b, 85.75.-d 



I. INTRODUCTION 

A. Electric generation of spin polarization 

Two purely electric mechanisms have been proposed to 
generate spin accumulations in electronic systems with 
spin-orbit interaction (SOI). They are current- induced 
spin polarization [l|, 0] and the spin Hall effect d, 
Both phenomena appear when a spin-orbit coupled con- 
ductor is traversed by an electric current. Current- 
induced spin polarization manifests itself in a weak bulk 
polarization of the electronic spins in a direction deter- 
mined by the electric current and the SOI, while the spin 
Hall effect manifests itself in a spin accumulation at the 
lateral edges of the sample. Though the resulting polar- 
izations are usually rather weak, both effects have been 
demonstrated experimentally p-Q . They offer an elegant 
alternative to spin injection and manipulation with fer- 
romagnets, Zeeman fields or by optical means, and pave 
the way toward spintronic devices with purely electric 
control of electronic spins. 

A dual aspect of these magnetoelectric effects is that 
they are accompanied by spin currents, flowing either 
inside the sample or leaving it. Early analytical investi- 
gations of these spin currents have been devoted to cal- 
culations of bulk spin conductivities, but it was soon re- 
alized that the latter are not always related to physically 
relevant quantities. Local approaches based on classi- 
cal kinetic equations have been used to investigate the 
conversion of magnetoelectric spin accumulations at the 
edge of a spin-orbit coupled sample into spin currents 
leaving the sample [1, [lol - [l^ . These quasiclassical the- 



ories have been quite successful, however they neglect 
coherent effects and, in their present form, arc not appro- 
priate to investigate local fluctuations and spatial corre- 
lations of spin polarization. Coherent mesoscopic effects 
on spin currents have recently attracted quite some the- 
oretical attention, both analytically jl5H20l | and numeri- 
cally , in both diffusive and ballistic systems, while 
so far only mesoscopic fluctuations of sample-integrated 
spin polarizations in the diffusive regime have been inves- 
tigated At present, little is known about local coher- 
ent fluctuations of spin polarizations, the scale on which 
they occur, how they are correlated, and how they in- 
fluence spin currents and their fluctuations, neither have 
current-induced spin polarizations in chaotic ballistic sys- 
tems been investigated so far. It is the purpose of this ar- 
ticle to start filling these gaps, by investigating the mag- 
netoelectric generation of spin polarizations in ballistic 
mesoscopic samples. 

To that end we construct a spin-probe model, which 
connects local spin-dependent chemical potentials to 
those imposed on a nearby probe. In this way, we are 
able to express local current-induced spin polarizations in 
terms of transmission coefficients deriving from the scat- 
tering matrix. Focusing on chaotic ballistic mesoscopic 
systems, we then use random matrix theory [25j to 
calculate the statistics of these polarizations. We find 
that, while the ensemble average polarization vanishes, 
it fluctuates from position to position (in a given sam- 
ple) and from sample to sample. For sufficiently large 
systems, the amplitude of these fluctuations by far ex- 
ceeds the magnitude of all current-induced polarizations 
investigated to date. Simultaneously, we show that this 
polarization is mostly uncorrelated with spin currents in 
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Figure 1: The spin-probe set-up. An open quantum dot with 
spin-orbit interaction is traversed by an electric current (/). 
The magnetoelectrically generated spin accumulation is mea- 
sured locally by an external probe (P) where spin-orbit inter- 
action is absent. The charge and spin-resolved electrochemi- 
cal potentials of the probe are set such that no spin nor charge 
current flow between the dot and the probe. When the probe 
is contacted to the system via a tunnel barrier of low trans- 
parency, we argue that the spin-dependent potentials in the 
probe can be identifled with the local spin-resolved potentials 
in the sample at the contact with the probe. This allows 
to relate spin accumulations to spin-dependent transmission 
coefficients. 

external leads and thus cannot easily be converted into 
a large spin current, at least not in linear response. 

B. The spin-probe model 

Fig. [T] gives a sketch of the system we focus on. We 
consider a two-dimensional chaotic ballistic quantum dot 
connected via reflectionless transport terminals to two 
external electron reservoirs. The reservoirs are biased 
and an electric current flows through the dot. Assum- 
ing that SOI is strong enough that the spin-orbit time is 
shorter than the mean electronic dwell time in the dot, 
Tso <^ td, spin rotational symmetry is totally broken in- 
side the dot. Together with chaotic boundary scattering, 
this leads to spin relaxation. Does the current flowing 
through this SOI coupled, spin-relaxing system generate 
a spin polarization ? On one hand, the relative positions 
of the contacts to the external electrodes determine an 
average electronic drift. In that sense, the situation is 
very similar to that considered by Edclstcin 0, and one 
would accordingly argue that a current-induced spin po- 
larization should emerge. On the other hand, however, 
electronic distribution functions are not well defined lo- 
cally in such a system where momentum and energy relax 
on the scale of the system itself. One might thus wonder 
how local spin-dependent potentials can be consistently 
defined. 

Here, we circumvent this difficulty by introducing an 
additional mobile terminal whose purpose it is to probe 



the local chemical potentials on the dot. We assume that 
spin currents are well defined in the probe as well as in 
the transport terminals, which presupposes that spin re- 
laxation processes are frozen there. In particular, SOI is 
present only in the dot. The probe is assumed to carry a 
single-transport channel. It is connected at one end to a 
macroscopic external reservoir, which defines its electro- 
chemical and spin-dependent chemical potentials unam- 
biguously, while its other end is contacted to the quan- 
tum dot via a tunnel barrier of transparency < Fp < 1 . 
When Fp <C 1, the probe does not perturb the dot and 
we can identify the spin-dependent chemical potentials 
on the dot with the corresponding probe potentials de- 
fined by the condition that neither char ge n or spin cur- 
rents fiow between the probe and the dot [26| . The probe 
being mobile, it can scan the whole dot area, and being 
single-channeled, it has maximal resolution. It can thus 
map the spin-polarization in the dot with maximal accu- 
racy. 



C. Summary of main results 

The spin-probe model becomes a powerful theoretical 
tool to calculate spin accumulations once electric and 
spin currents in the probe and the two transport termi- 
nals are related to chemical potentials via the scattering 
approach to transport [23, ll^. When the probe car- 
ries a single transport channel, time-reversal symmetry 
and the unitarity of the scattering matrix ensure that 
the reflection coefficients from the probe back to itself 
are diagonal in spin indices, independently of Fp. This 
allows to directly express the electrochemical and spin- 
dependent chemical potentials of the probe in terms of 
spin-dependent transmission coefficients from the trans- 
port leads to the probe, and the chemical potentials at 
those leads. This is done below in Eq. We addition- 
ally find that, for a single-channel probe, the dependence 
that transmission coefficients have on Fp factorizes. This 
has the important, and somewhat surprising consequence 
that the chemical potentials ensuring the vanishing of the 
currents through the probe, being given by the ratio of 
transmission coefficients, do not depend on Fp. This is 
an exact result, valid as long as one has a single-channel 
probe, and we stress that it is true not only on average 
nor to some leading order. Consequently, we are able to 
obtain the value of the spin-dependent potentials on the 
probe in the limit Fp ^ 0, where it is justified to identify 
the potentials on the dot with those on the probe, while 
performing the calculation for the analytically simpler 
case Fp = 1. 

Our scattering theory of spin accumulations can be ap- 
plied to both diffusive and ballistic mesoscopic systems 
in the linear response regime. Because current-induced 
spin polarization has not been considered in these sys- 
tems so far, we specialize here to ballistic systems, and 
accordingly apply random matrix theory (RMT) [Ij, HH 
to calculate average and fiuctuations of the spin-resolved 
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local chemical potentials. We perform our calculation 
for Tp = 1 to avoid unnecessary complications arising 
from the presence of a Poisson kernel for Fp < 1 [23 |. 
We find that, within RMT, there is no current-induced 
local spin polarization on mesoscopic average. However, 
sample-dependent or probe position-dependent fluctua- 
tions make local spin-dependent chemical potentials fluc- 
tuate with a variance 
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given by the voltage bias V between the two transport 
leads and the number N of channels each of them carries. 
The above formula gives the leading order contribution in 
1/N 1. The exact formula is given below in Eq. ([T4|. 
Within RMT, the variance is the same for (3 — x,y and 
z. 

The transmission coefficients from the transport ter- 
minals to the probe are elements of a scattering matrix 
which depends on the position of the probe. Therefore, 
as the probe scans the sample, these coefficients fluctu- 
ate. Below, we conjecture and numerically confirm that, 
within RMT, moving a single-channel terminal by a sin- 
gle Fermi wavelength is sufficient to decorrelatc transmis- 
sion coefficients to and from that terminal. This has for 
consequence that the total polarization integrated over 
the entire dot, while still averaging to zero, acquires a 
sample-to-sample variance given by the single-channel lo- 
cal variance multiplied by the total number kpA/Air of 
uncorrelated probe positions. Thus, 
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fraction of fcp. It is also much larger than the typical 
Duckheim-Loss polarization [2^, 
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predicted for kaL ^ 1. Despite its magnitude, (S'jyp) 
cannot easily be converted into a large spin current. Be- 
low we demonstrate this by calculating the correlator of 
the spin polarization with the spin current in external 
leads, and show that this correlator vanishes unless the 
probe is within a Fermi wavelength of the contact to the 
terminal. In other words, RMT predicts that only ac- 
cumulations right at the boundary of the sample can be 
extracted and converted into spin currents in the linear 
regime, which is consistant with the universal fluctua- 
tions of spin currents found in Refs. [l6l - [l9| . It is at 
present unclear if more spin current can be extracted in 
nonlinear regimes of transport. 

Having summarized our method and main findings, 
we now proceed to present our theory more formally. 
In the next Section [Hi we construct a scattering the- 
ory of current-induced spin polarization. In Section IIIIl 
we focus on ballistic systems and apply RMT to calcu- 
late the mean and variance of the local as well as of the 
sample-integrated polarizations. We further demonstrate 
the power of our approach by calculating correlators of 
spin polarizations with spin currents. Section IIVI is de- 
voted to a numerical proof-of-principle of our theory and 
a quantitative conflrmation of our main result, Eq. ([3]). 
We conclude with a brief discussion of our results and of 
future challenges and demonstrate the Fp-independcnce 
of the probe voltages for Np = 1 in an appendix. 



Because fi[^l gives the chemical potential difference be- 
tween the two possible spin orientation along the spin 
axis /3 = x,y,z, dividing it by twice the level spacing 
converts it into a number of polarized spins in that di- 
rection. This gives the typical number of polarized spins 
in a given sample as 
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with the dot level spacing 6 = 2n/mA^ the width W of 
the transport terminals and the linear dot size L. Eq. ([3]) 
is one of our main results. It predicts that linear electric 
transport across a typical, ballistic mesoscopic quantum 
dot generates a giant polarization of the electronic spins 
inside the dot in the semiclassical limit k-pL ^ 1, and 
for a large enough dwell time L/W ^ 1. This typical 
polarization is in particular much larger than the average 
magnetoelectric polarization predicted by Edelstein Q 
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in a diffusive dot with clastic mean free path £ <^ L and 
with Rashba SOI whose strength ka is usually only a 



II. SCATTERING THEORY OF 
CURRENT-INDUCED SPIN POLARIZATION 

In its standard form, the scattering approach to trans- 
port gives a linear relation between currents and voltages 
measured in external leads [l^, HI] . It applies equally 
well to diffusive and ballistic systems. Extended to take 
account of spin currents and chemical potentials in a mul- 
titerminal geometry, the relation reads 
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where Greek letters label spin indices and Latin let- 
ters label electric terminals. Here, Ni gives the num- 
ber of transport channels in terminal z, F^ the trans- 
parency of the contact between that terminal and the 
dot, ^ = + eV^ gives the electrochemical potential 
in terms of the Fermi energy E-p and the applied volt- 
age Vi, and give the chemical potential difference 
between the two possible spin orientation along the spin 
axis P = x,y,z [2§|. We further defined = /; and 
as the charge current and the three components of 
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the spin current vector respectively, and introduced the 
spin-dependent transmission coefficients [l9j 

^M)^ E Tr[4„a-Wa^], (7) 

where ct", a = x,y,z are Pauli matrices and cr^ is the 
identity matrix. The trace in Eq. ([7]) is taken over the 
spin degree of freedom and is a 2x2 matrix of spin- 
dependent transmission amplitudes from channel n in 
lead J to channel m in lead i. 

We are interested in the mesoscopic counterpart of 
the magnetoelectric geometry of Ref. 0], and consider 
a quantum dot connected to two external transport ter- 
minals L and R. An electric bias voltage V is applied 
between the terminals, /Li^°^ — = ^A^h^ + Ey — 6^/2, 
where we further assume that there is no spin accumu- 
lation, = 0, /3 ^ 0. To measure local chemical 
potentials, we introduce an external probe with control- 
lable chemical potentials, which scans the dot surface. 
The situation is sketched in Fig. [T] The chemical poten- 
tials in the probe are determined by the condition that 
no charge nor spin current flows through it, i.e. they are 
solutions of 

Ip^-l E E(2rpA^^'5.,<5p,-T(ff);.f ^0, 

j=L.R,P P 

(8) 

and depend a priori on the transparency of the contact 
between the probe and the dot. Eq. ([8|) determines the 
local chemical potentials ^jl^P as a function of the elec- 
trochemical potentials /i^°p imposed on the longitudi- 
nal leads. We consider a single-channel probe A^p = 1, 
which not only delivers the sharpest resolution of the 
dot's chemical potentials landscape, but moreover en- 
sures that the spin reflection matrix T^p^ from the probe 
to itself is exactly (i.e. not only on average) diagonal in 
spin indices. This follows from time-reversal symmetry, 
current conservation and gauge invariance [30| , which ad- 
ditionally require T^p-* = Tpp , Va. One then straight- 
forwardly obtains the chemical potentials in the probe 
as 

,f - - «]eV2[<^ + (9) 

In particular, for Nl,r ^ 1, the electrochemical potential 
in the probe is given by 

Mli'^ = [T<°?-T&l^/4Fp, (10) 

which reproduces a result first derived in Ref. [3l|. In 
diffusive samples, this gives a linear average decay of /Xp"* 
from eV/2 to —eV/2 as the probe is moved across the 
sample from regions more easily accessible from the left 
lead (where > ) to regions which are more 

easily accessible from the right. 

We next make the key observation that the Fp- 
dependence of the transmission coefficients T^'^^ between 



the probe and the transport terminals I = L,R exactly 
factorizes when the probe is single-channeled. This re- 
markable property is apparently mentioned here for the 
first time and is discussed in some detail in Appendix A. 
Consequently, the chemical potentials in Eq. ([9]) remain 
the same, regardless of the transparency of the tunnel 
barrier between the dot and the probe, when the lat- 
ter carries only a single channel. Our aim being to ob- 
tain the current-induced spin polarization and thus the 
spin-dependent chemical potentials on the dot, we are 
quite naturally interested in the limit F — >■ 0, where the 
probe does not perturb the dot and accordingly the spin- 
dependent chemical potentials on the dot equal those on 
the probe. The observation we just made allows us to 
access that limit via an analytically simpler calculation 
at Fp = 1 which we present in the next section. 

To get the total spin polarization, we scan the probe 
over the whole sample. The transmission amplitudes tmn 
in Eq. ([7]) are elements of a scattering matrix that de- 
pends parametrically on the probe position. For a wide 
probe, Np ^ 1, we argue semiclassically that the corre- 
lator (T^f ^T^°i^) - (Tj?/')(Tj."^') decays to zero when 
the distance between P and P' is of order the width of 
the probe. This is so, because beyond that, correlated 
escapes of pairs of trajectories no longer contribute [13 ■ 
We next make the leap of faith that this semiclassical 
argument can be extended to the case of present inter- 
est of a single-channel probe, for which we accordingly 
conjecture that two measurements performed a distance 
^ kp^ apart deliver uncorrclated outcomes. To the best 
of our knowledge, this problem has not been investigated 
yet, and we will numerically confirm the validity of this 
conjecture below. Further assuming a sharp decay of the 
correlator at k^^, we obtain that the sample- integrated 

chemical potentials are given by the sum ^[^l = ^p /ip ^ 
over kpA/Air uncorrclated probe positions on the surface 
of our two-dimensional quantum dot. 



III. RANDOM MATRIX THEORY 

Having expressed local and sample-integrated polar- 
izations in terms of transmission coefficients, Eq. 0, we 
next proceed to evaluate them quantitatively. We focus 
on ballistic chaotic cavities, and accordingly use RMT 
to calculate the average and variance of fj.^ and /ijft. 
This calculation is easier for Fp = 1, since then the 
ensemble-average scattering matrix vanishes [23|. Be- 
cause of the factorization property of the F p-dcpendcnce 
of transmission coefficients mentioned above, the spin- 
dependent potentials on the dot arc equal to those on 
the probe for Fp — !■ 0, which in their turn remain un- 
changed as Fp — > 1. Thus we can extract the dot's 
spin-dependent potentials fix^J at point xp underneath 
the probe from a calculation at Fp = 1, because fii^J = 
fj,p\Tp — >• 0) = /Zp^(Fp — 1). Because we are inter- 
ested in systems with time reversal symmetry but totally 
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broken spin rotational symmetry, we take the scattering 
matrix as an element of the circular symplectic ensem- 
ble of random matrices [1^. RMT does not depend on 
the exact form of SOI, but only requires that SOI totally 
breaks spin rotational symmetry. Our results are thus ex- 
pected to apply for Rashba and Dresselhaus SOI equally 
well, as long as one has a ballistic chaotic quantum dot 
with totally broken spin rotational symmetry. 
For Np = 1, we get from Rcf. [l^ 



/RMT 



Nl + Nr + 1/2 



from which wo conclude, together with Eq. (|9]), that 



(11) 
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(Mtot)RMT 



0. 
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RMT thus predicts the vanishing of the spin polarization 
on average. The latter is understood in the standard 
mcsoscopic sense, as an average over changes in the dot 
shape, the overall chemical potential or the position of 
the leads. 

Chemical potentials nevertheless exhibit spatial fluctu- 
ations in a single sample, whose magnitude can be esti- 
mated by calculating var /Xp ■* . Of further interest are also 
mesoscopic fluctuations of the sample-integrated chemi- 
cal potentials whose magnitude wc analyze with var ^[^l . 

To calculate the variance of in Eq. we take the 
variance of the numerator and divide it by the squared 
average of the denominator, a procedure which neglects 
subdominant corr ections 0(7V^3^). The variance of the 
spin-dependent transmission coefficients have been cal- 
culated in Ref. Ilg. For i = L, R, wc have 



1) - NpNf 



covar(<°\Tlfp°)) = 



Nt{2Nt - 1){2Nt - 3) 
-ANpNlNr 



Nt{2Nt - l)(2iVT - 3) 



, (13a) 
, (13b) 



with Nt ^ Nl + 
we obtain 

(/3) / eV 

var /ip' = 16 



Np. For our case with Np = 1, 



NlNr 



7VT(2iVT - l)(2iVT - 3) 
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For = TVp = TV > 1, Eq. 1^ goes into Eq.©. 

With the further assumptions, discussed in the pre- 
vious section, that (i) spin relaxation processes are 



frozen in the probe and (ii) the correlator C^pl 7p,j^ 



/3)T-(a^')\ 



{T^Pi}){T^pPi!) sharply decays when the distance between 
the probe positions P and P' reaches the Fermi wave- 
length, we get the typical total spin-dependent chemical 
potentials 



polarization, we convert them into the variance of the 
spin polarization by dividing them with 25. For Nl = 
Np^ = N, wc obtain the variance of the total number of 
polarized spins on the dot as 



var S 
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(2iV+l)(4Ar + l)(4iV-l) 
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which, for — > oo, gives Eq. ([3]). The validity of 
Eqs. ([T414T51) is numerically confirmed below in Fig. [3] 

The power of our scattering theory is perhaps even 
more evident when one realizes that it allows to calculate 
correlation functions between local spin polarizations in- 
side the sample and spin currents in the leads. The spin 
current in lead L is given by 



(a) 



h 



(17) 



and each of these three terms gives a separate contri- 
bution to the correlator (/["Vp^) - (4"^)(Mp^)- For 
a = x,y, z, the contribution of the third term van- 
ishes because it contains the average of three spin- 
dependent transmission coefficients. Furthermore, RMT 
gives covar(T^"°\T]f°'') cx SikSapSao if j ^ P- Therefore, 
for a ~ x,y, z, the first two terms do not contribute ei- 
ther, unless the probe is placed in the immediate vicinity, 
i.e. within a single Fermi wavelength of the contact to 
the L lead. In this case, nonvanishing covariances appear 
in ((/{"•'/ip ■*)); see Eq. (jl3bp . and one obtains 



((/rV^«)) 



-SSq^^Sp^^NlNr e'V^ 
Nt{2Nt ~ 1){2Nt - 3) h 



(18) 



where Sp^^, requires that the probe be placed within fcp ^ 
of the boundary between the SOI-coupled dot and the 
lead. Summing over all possible probe positions, one 
obtains a universal correlator )) = 0{N'^) for 

Nl.r ~ iV 3> 1, which is consistent with the universality 
of spin currents fluctuations [l6l - [l9| . RMT thus predicts 
that coherent polarizations correlate with the spin cur- 
rents only if they are right at the sample exit. 



IV. NUMERICS 

We numerically illustrate the validity of our scattering 
theory of spin polarizations and confirm our main predic- 
tion, Eqs. (|14m5p . We use tight-binding models defined 



by the Hamiltonian 



(19) 
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var/iiot 



— — var Up ' , 
47r 
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As justified above, we identify these potentials with those 
on the dot. Assuming that they arc generated by a spin 



where the operators cj^ (qs) create (destroy) a spin-1/2 
fermion with spin s =t,i and (t"^, indicate elements of 
Pauli matrices. To give a proof-of-principlc for our ap- 
proach, we consider a diffusive model, while quantitative 
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Figure 2: Numerical proof-of-principle for the local measurement of spin accumulations on a single, diffusive, Rashba-coupled 
lattice of size 20 x 50 traversed by an electric current in the positive a;-direction (from left to right). The on-site disorder 
and Rashba spin-orbit interaction correspond to a mean free path £ = 6a and a spin-orbit length /so = 9a. The system has 
hard-wall boundary conditions at its top and bottom boundaries. From left to right, spin accumulations /i'''' with jS = x, y 
and z. Only displays a finite average magnetoelectric effect (central panel), as predicted in Ref. 0]. Fluctuations exist for 
all spin orientations. The color scale is in units of the voltage bias eV . Note the different color scales in different panels. 



confirmations of our tlicory are solely based on the RMT 
model. 

For the diffusive model, i = {ixjiy) labels a site on 
a two-dimensional lattice of size Lx Ly ~ {Nxo) x 
[Nyo) with spacing a, giving the size M ~ Nx x Ny 

of the Hamiltonian matrix. We further take /i,-"' = 

itsoSa,xSi^J^ {Siy+a,jy + Si^—a.jy ) — itso5a,y5iy,jy {^i^+aj^ + 
Sij:-a,jj;)~'tSa,o[Si^,j^ {diy+a,jy ~\~Si^-a,jy) + 5i^jy (Si^+aJ^ + 

Si^-a.j^)] + WiSa.oSij, corresponding to an on-site 
disorder potential with randomly distributed Wi S 
[— u'/2, w/2], a nearest-neighbor hopping t and a Rashba 
spin-orbit interaction tgo- These parameters give a spin- 
orbit length ^so ~ TToi/iso and an elastic mean free path 
£ « iSat^/w^ close to half filling. 

For the RMT model, we takeY^^fa"' to be 
randomly distributed quaternions [24 . |25| . with time- 
reversal symmetry requiring hf^'^ = h^^^ , and /il"-* = 

—h^ji"^ for a = x,y,z. Generally, the indices i,j have 
no particular physical meaning in RMT and the size M 
of the Hamiltonian matrix is a free parameter. Here we 
attach the system to external leads, and also consider i 
and j to label sites on a two-dimensional lattice whose 
Hamiltonian has an infinite-range random hopping. Once 
the sharp decay of spin-spin correlations is confirmed, 
M = A:p^/47r gives the total number of uncorrelated 
probe positions. For both models, all our simulations 
fix the Fermi energy slightly below half-filling. 

To construct the scattering matrix, we follow the 
method of Ref. [11]. We introduce a M x Nt projection 
matrix W which connects the system to external leads. 
This matrix has only one nonzero element per row, at 
a spot determined by the coupling between the system 
and the leads. For an ideal contact, all nonzero matrix 
elements of W are unity. The A^t x scattering matrix 
§ associated with J{ is given by 

^^l-2'KiW'^{EF -'K + iTiWW'^y^W , (20) 

with the overall Fermi energy E^. The transmission co- 



efficients of Eq. ([7]) are straightforwardly obtained from 
S. For the diffusive model, this method is certainly not 
the most efficient numerically. We still rely on it, because 
for the RMT model, the infinite-range hopping renders 
the implementation of recursive Green's function algo- 
rithms problematic. Note that it is necessary to first 
generate random Hamiltonians from which random scat- 
tering matrices are constructed in order to consistently 
define probe positions and investigate whether spin-spin 
correlations decay sharply beyond a Fermi wavelength. 
To investigate the decay of these correlations, we gener- 
ate a random "K and take 

(In, \ 

W = \ VM-N,-Nn , (21) 

V J 

where for each choice of W, only a single component of 
the column vector v (whose index indicates its length) is 
nonzero and equal to one. The position of that nonzero 
component determines the position of the probe on the 
two-dimensional lattice. Our Hamiltonians are defined 
on rectangular two-dimensional lattices, and we contact 
the transport modes at two of its edges, which is repre- 
sented by the two identity matrices in Eq. (HI]). 

We first show numerical proof-of-principle for the spin- 
probe mode. In Fig. [2J we present data for the local 

spin polarizations /Xp'', /? = x,y,z. The three panels 
clearly show that, for transport in the a;-direction and 
a Rashba SOI, there is a finite average magnetoelectri- 
cally induced polarization along the y-axis only, as pre- 
dicted in Ref. 0- Furthermore, there are local (as well 
as sample-to-sample) fluctuations of the spin polarization 
in all directions, as predicted in Ref. j23| . All accumu- 
lations vanish on average in the RMT model, moreover 
spatial fluctuations occur on the smallest possible scale 
of a single site, as expected, and we therefore do not show 
them here. 

We next conflrm our conjecture of a sharp decay of 
correlations, and to that end, we calculate the spin-spin 
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M N /(2N+1)(4N+1)(4N-1) 



Figure 3: Main panel: numerically obtained variance of the 

(x) 

sample-integrated spin-dependent potential nl^l in the RMT 
model (circles) vs. the theoretical prediction of Eq. p6p 
(dashed line) for Nl = Nr = N e [3,25] and M G [64,625]. 
Deviations from the theory occur at smaller number of chan- 
nels, for which fluctuations in the denominator of Eq. (|9]) can 
no longer be neglected. The behavior of fi[^l and is the 
same as the one of Mtot- Inset; spin-spin correlation func- 
tion D{n) = {S^'^\i)S^^\i + na)) for the RMT model (black 
circles) and the diffusive model (red squares) with £ = 8a, 
lao = 12a. In both cases, we find that the correlator does not 
depend on the spin orientation /3. 



V. CONCLUSIONS 

Wc have demonstrated how introducing an external 
probe allows to relate local current-induced spin polar- 
izations in a mesoscopic system to spin-resolved transmis- 
sion coefficients. Focusing on ballistic systems, we used 
random matrix theory to compute the statistics of these 
polarizations. We predicted that a ballistic quantum dot 
with large contacts to external transport leads - but still 
large dwell time - typically carries a spin polarization 
that exceeds those previously investigated by a paramet- 
rically large factor [see Eqs. ([SHS])]. Perhaps most impor- 
tantly, we used this spin-probe model to calculate for the 
first time quantum correlations between local spin polar- 
izations and spin currents, beyond the classical Fick law. 
We showed that the reported large spin polarizations cor- 
relate only weakly with spin currents, thus they cannot 
be easily converted into spin currents - and vice-versa. 
As disappointing as this result may be, it is consistent 
with the universality of spin current fluctuations [l6l - [l9| . 

The reported accumulations are coherent in nature, 
and one might wonder how they disappear as the system 
is coupled to external sources of noise. In this respect 
it is interesting to note that for a multichannel probe, 
Eq. ([HI) becomes, for N,Np » 1, 



var Si 



(/?) 



kpA 
UAttN 



eV 
'26 



1 



(1 + Np/2Nf 



(22) 



correlation function D{n) ~ {S^^\i)S^^\i + nae)). For 
both diffusive as well as RMT models we found no sig- 
nificant anisotropy, i.e. D was the same for e in the x 
or y direction. This is of course expected for the RMT 
model, but would deserve some more investigations in 
the diffusive case. The Fermi wavelength being equal 
to the lattice spacing close to half filling, one expects 
D{n) = (5„o for the RMT model. This is confirmed 
in the inset of Fig. [31 where we additionally see that 
D{n) is also sharply peaked at n = in the diffusive 
model, with some correlation persisting over few, ^2-3 
sites. This sharp decay of D{n) in the diffusive model 
suggests that our conjecture of a fast decay of the corre- 
lator (T^f^T^"^^) - (T^t'^)(^p"L ) applies beyond RMT. 
This also deserves further investigations. 

Finally focusing on the RMT model. Fig. [3] confirms 
our main results, Eq. (jl4j|15|) . by showing the scaling 
of the variance of the sample-integrated spin-dependent 
chemical potential in symmetric samples Nl = Np = N 
as a function of N and the size M of the random matrix 
Hamiltonian. Small deviations from our analytical pre- 
dictions occur at small N for which our procedure of re- 
placing the denominator in Eq. ^ by its average squared 

when calculating var /ip ^ is no longer justified. All these 
results numerically confirm our scattering theory of spin 
polarizations. 



Enlarging the spin-probe damps the variance of the polar- 
izations in a similar wayas a dephasing probe damps con- 
ductance fluctuations |34l | , with however a different para- 
metric dependence on Np/N. Note that the probe gen- 
crating this damping is spin-conserving on time-average, 
as there is no spin current flowing through it. 

Beside decoherence effects, there are a number of in- 
teresting open problems, related for instance to spatial 
fluctuations of spin polarizations and our assumption of 
a sharp decay of correlations at fcp^. Finding out how 
the polarizations correlate to spin currents in diffusive 
systems, as well as in true ballistic systems treated be- 
yond RMT is also of interest. Ways to extract the large 
typical polarizations we reported in the nonlinear trans- 
port regime or by optical means as in Ref. Q should also 
be explored. 
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Appendix A 

We demonstrate that the potentials ensuring zero cur- 
rent through a single- channel probe do not depend on 
the transparency of the contact between the probe and 
the system. According to Eq. ([9]), the chemical poten- 
tials in an external probe connected to a voltage-biased 
two-terminal cavity are given by the ratio of the differ- 
ence and the sum of transmission coefficients between 
the transport leads and the probe. We show that the 
Fp-dependence of the coefficients 7''pf^ factorizes when 
Np = f so that Fp docs not affect the chemical potential. 

I 



Our starting point is Eq. (|20|) for the scattering matrix, 
which wc rewrite here, 

§^l-2T:iW\EF -^ + iTrWW^)-^W . (23) 

The matrix WW^ = VV^ + PP^ projects onto the sys- 
tem's states coupled to the transport leads (T^T^^) and 
the probe {PP"^). We can write a perturbative series in 
PP"^ for the elements of the scattering matrix connect- 
ing the leads to the probe {ji is a channel index on the 
I = L,R leads) 



Spj, = -2Tri"^P*pV^>j,{x\{EF -3< + i7TWW'')-^\x') 
= -2TTi^P,*pV^>j, [{x\{Ef - 0< + iTTVV'')-^\x') 

x,x' 

-iT:{x\{EF-:K + inVV^y^PP''{EF-:K + i7:VV^y^\x) + ...] . (24) 

The transparency of the contact between probe and system is contained exclusively in the matrix P. We are free to 
choose any orthogonal basis {\x)} for the system and pick one where the M x 1 coupling matrix between probe and 
system is Pxp ~ uopSxP- One then has fo r the M x M matrix {x\PP'^\x') = {P\PP''\P)SxpSx'p = \ljp\ SxpSx'p, with 
|wp|2 = {]\I5/n^Tp){2 - Fp - 2Vf -Fp) [H. One gets 

OO 

Spj, = -2niuj;^{-nT\ujpf {P\{Ef -'K + inVV^)-^\P))" xJ^'^^'jAPHEf - ^ + ittVV^)-^\x') 

n=0 x' 

= f{rp)xJ2Vx'jAP\{EF-^ + i7rVV^r'\x'). (25) 



One sees that, because the probe is single-channeled, 
the Fp-dependcnt part of §pjj factorizes, and the same 

obviously holds for transmission coefficients Tpj'' 

l§PjjP- This is not the case for a multichannel 
probe, in which case \P) is no longer unique. Thus, 
for Np = 1, the electrochemical potential in Eq. (fTO|) 



do not depend on the strength of the tunnel barrier. 
The argument is straightforwardly extended to the spin- 
dependent chemical potentials in Eq. (j9l), because for a 
single-channel probe, the reflection coefficient T^p is di- 
agonal in spin-indices. 



[1] F.T. Vasko and N.A. Prima, Sov. Phys. Solid State 21, 
994 (1979); L.S. Levitov, Y.V. Nazarov, and G.M. Eliash- 
berg, Sov. Phys. JETP 61, 133 (1985); A.G. Aronov and 
Yu.B. Lyanda-Geller, JETP Lett. 50, 431 (1989). 

[2] V. M. Edelstein, Solid State Comm. 73, 233 (1990). 

[3] M.I. Dyakonov and V.I. Perel, Sov. Phys. JETP Lett. 13, 
467 (1971). 

[4] For a review on the spin Hall effect, see: H.-A. Engel, 
E.I. Rashba, and B.I. Halperin, Theory of Spin Hall Ef- 
fects in Semiconductors, Handbook of Magnetism and 
Advanced Magnetic Materials, Wiley and Sons (2007); 
,arXiv:cond -mat / 0603306 

[5] Y.K. Kato, R.C. Myers, A. C. Gossard, and D. D. 



Awschalom, Science 306, 1910 (2004). 
[6] J. Wunderlich, B. Kastner, J. Sinova, and T. Jungwirth, 

Phys. Rev. Lett. 94, 047204 (2005). 
[7] S.O. Valenzuela and M. Tinkham, Nature 442, 176 

(2006). 

[8] E. Saitoh, M. Ueda, H. Miyajima, and G. Tatara, Appl. 

Phys. Lett. 88, 182509 (2006); T. Kimura, Y. Otani, T. 

Sato, S. Takahashi, and S. Maekawa, Phys. Rev. Lett. 

98, 156601 (2007); 98, 249901(E) (2007). 
[9] H. Zhao, E.J. Loren, H.M. van Driel, and A.L. Smirl, 

Phys. Rev. Lett. 96, 246601 (2006). 
[10] E.G. Mishchenko, A.V. Shytov, and B.I. Halperin, Phys. 

Rev. Lett. 93, 226602 (2004). 



9 



[11] A.A. Burkov, A.S. Nunez, and A.H. MacDonald, Phys. 

Rev. B 70 155308 (2004). 
[12] A.G. Mal'shukov, L.Y. Wang, C.S. Chu, and K.A. Chao, 

Phys. Rev. Lett. 95, 146601 (2005). 
[13] i. Adagideli, M. Scheid, M. Wimmer, G.E.W. Bauer, and 

K. Richter, New J. Phys. 9, 382 (2007). 
[14] Y. Tserkovnyak, B.I. Halperin, A.A. Kovalev, and A. 

Brataas, Phys. Rev. B 76, 085319 (2007). 
[15] O. Chalaev and D. Loss, Phys. Rev. B 71, 245318 (2005). 
[16] J. H. Bardarson, I. Adagideh, and Ph. Jacquod, Phys. 

Rev. Lett. 98, 196601 (2007). 
[17] Y. V. Nazarov, New J. Phys. 9, 352 (2007). 
[18] J. J. Krich and B. L Halperin, Phys. Rev. B 78, 035338 

(2008). 

[19] I. Adagideli, J.H. Bardarson, and Ph. Jacquod, J. Phys.; 

Condens. Matter 21, 155503 (2009). 
[20] J.J. Krich, Phys. Rev. B 80, 245313 (2009). 
[21] B.K. Nikolic, L.P. Zarbo, and S. Souma, Phys. Rev. B 

72, 075361 (2005). 
[22] W. Ren, Z. Qiao, J. Wang, Q. Sun, and H. Guo, Phys. 

Rev. Lett. 97, 066603 (2006). 
[23] M. Duckheim and D. Loss, Phys. Rev. Lett. 101, 226602 

(2008). 

[24] P.W. Brouwer and C.W.J. Beenakker, J. Math. Phys. 37, 
4904 (1996). 

[25] M.L. Mehta, Random Matrices, Academic Press, London 
(1990). 

[26] Because of Coulomb interactions, such a direct relation 
does not necessarily exist for the electrochemical poten- 
tial on the dot, see e.g.: M. Biittiker, J. Phys.: Condens. 
Matter 5, 9361 (1993). 

[27] M. Biittiker, Phys. Rev. Lett. 57, 1761 (1986). 



[28] Y. Imry in Directions in Condensed Matter Physics, G. 
Grinstein and G. Mazenko eds., World Scientific, Singa- 
pore (1986). 

[29] For /3 = x,y, z, y!f^ gives the chemical potential differ- 
ence between the two possible spin orientations along the 
corresponding spin axis, while = Ep + eVi gives the 
usual electrochemical potential. We refer to electrochem- 
ical and spin-dependent potentials collectively as "chem- 
ical potentials". Note the plural. "Chemical potential" 
still retains its standard meaning. 

[30] F. Zhai and H.Q. Xu, Phys. Rev. Lett. 94, 246601 (2005); 
A.A. Kiselev and K.W. Kim, Phys. Rev. B 71, 153315 
(2005). 

[31] M. Biittiker, Phys. Rev. B 40, 3409 (1989). 

[32] Correlated escape of pairs of classical trajectories and its 
infiuence on semiclassical transport was first discussed 
in: Ph. Jacquod and R.S. Whitney, Phys. Rev. B. 73 , 
195115 (2006). The only nonvanishing semiclassical con- 

tributions to {^'t;'^^^^'^) - for a ^ have 

trajectories with correlated escape; Ph. Jacquod, unpub- 
lished (2008). 

[33] C. Mahaux and H.A. Weidenmiiller, Shell-Model Ap- 
proach to Nuclear Reactions, North-Holland, Amsterdam 
(1969). 

[34] M. Biittiker, Phys. Rev. B 33, 3020 (1986); H.U. 
Baranger and P.A. Mello, Waves in Random Media 9, 
105 (1999). 

[35] Th e eigen value of PPt |^^|2 ^ [M5/-k^Vp){2 - Fp ± 
2-^1 — Fp) is not uniquely determined by Fp. Here we 
take the smallest of the two solutions, for which the ex- 
pansion of Eq. (|24[) is guaranteed to converge. 



